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INTRODUCTION 


The purpose of this study was to select monthly weather variables that 
could be used to estimate wheat yields for the wheat growing areas of Brazil. 
Wheat is grown in seven states in Brazil: Rio Grande do Sul, Santa Catarina, 
Parana, Sao Paulo, Mato Grosso, Goias, and Minas Gerais. Rio Grande do Sul 
is Brazil's original wheat-producing state; until 1972 it was the country's 
most important production area. Increasingly, however, Parana, Mato Grosso 
and Sao Paulo have grown in importance. In 1977, Parana took over the number 
one spot. Figure 1 shows the wheat-growing areas of Brazil. 

Although wheat has been grown since the sixteenth century, Brazil has yet 
to develop a high-quality wheat variety that produces well under the country's 
widely varied climatic conditions. Brazil's wheat crops have continuously 
been plagued with problems resulting in consistently low yields. Frosts and 
plant disease occasionally reduce yields. The high cost discourages use of 
fungicides. Expansion of acreage sown to wheat was met with cultivating 
problems and high costs. New land areas cultivated in wheat are highly acidic 
and low in fertility. Fertilizer is costly. Late-season rains frequently 
delay harvest and reduce yield. 

The climate of the southern wheat-producing states is "subtropical humid" 
rainfall is relatively abundant and well-distributed throughout the year, with 
slightly more rainfall in the warm months. There is usually no season of 
drought. The northern states are "semi arid" with a winter dry season and less 
total annual rainfall. Summers are hot and winters are mild. Parana and 
southern Sao Paulo are the northern limit for frost ocurrence. 

Wheat is planted in the months of April through June and is harvested in 


November and December 
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Figure 1. Wheat-growing areas of Brazil 

(J. McQuigg, R. Willis, 1982, 'personal 

communication) 
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METHOD 


Multiple regression analysis of yield with selected agroclimatic indices 
was used to derive a suitable model. The index P-PET (precipitation minus 
potential evapotranspiration) was used in the regression equations to represent 
available soil moisture, monthly precipitation, and monthly maximum temperature. 

The regression equation is: 

A 

Y = CX + BjTX-j + B2Ri + B3 (P-PET)j + E 


where 


A 
Y = 

Estimated yield. 



a = 

Constant, 



B j * 

Coefficients of variables j = 

1-3, 

TXi = 

Maximum temperature for 

month 

i. 

Ri = 

Total precipitation for 

month 

i. 

(P-PET)-j = 

Precipitation minus PET 

for month i 

E = 

Unexplained error. 




In developing the model, various procedures of the Statistical Analysis 
System (SAS Institute, Inc., 1979) were used. The procedures used and the 
operations performed with each are summarized in the Appendix. The selected 
model had the highest r 2 and included variables that were significant at the 
10 per cent level and agronomical ly meaningful. 

DATA 

The Brazil crop data were obtained from the Foreign Agricultural Service 
(Sam Ruff, Personal Communication, 1982). The data was recorded with year 
of yield as year of harvest, so the weather influencing the crop occurred 
during year - 1. 
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Meteorological data from 1972 through 197? were used to model because 
there Is no apparent trend In the yield data during this period. Furthermore, 
1977 represented the latest available data. Table 1 lists the stations used to 
derive the meteorological data sets for each state. Figure 2 shows the 
location of each station. 


PROCEDURES 

Since Rio Grande do Sul's meteorological data has the longest period of 
record. Initial models were developed for this state alone. Various weather 
variables were tried in regression equations In many different combinations. 

The results were not good. Several different trend, were tried with different 
regression models, but none were significant at the 10 per cent level. From 
historical accounts, it Is believed that some sort of trend of Increased yield 
began In 1962 because of Increased use of fertilizer and more adaptable 
varieties of wheat. Yet the yield data did not indicate this. It was decided 
to model for years 1972-1979 for which it was believed no trend existed. 
Eliminating years of data created the problem of fewer degrees of freedom. 
However, by combining data for the states of Rio Grande do Sul, Parana, and 
Santa Catarina, a covariance model could be devloped. 

Variables used in the regression equations for the covariance model. 

Included "dummy variables" for Parana and Santa Catarina. The "dummy variables" 
adjust the contributions to yield of both states to a base yield which. In this 
case. Is Rio Grande do Sul's yield. The "dummy variable" for Parana was not 
significant at the 10 per cent level in the final model. The coefficient for 
Santa Catarina's "dummy variable" is negative, indicating that its yield is 
below that of the norm set by Rio Grande do Sul. 
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STATE 

METEOROLOGICAL STATION 

WMO NUMBER 

Mato Grosso 

Campo Grande 

83611 


Ponta Pora 

83702 


Pres Prudente 

83716 

Parana 

Londrina 

83766 


Santa Branca 

83810 


Curitiba 

83842 


Porto Uniao 

83864 

Rio Grande do Sul 

Ijui 

83805 


Sao Borja 

83901 


Passo Fundo 

83914 


Julio de Castilhas 

83935 


Santa Maria 

83936 


Sao Gabriel 

83957 


Bage 

83980 

Santa Catarina 

Porto Uniao 

83864 


Vacaria 

83918 


Sao Joaquin 

83920 

Sao Paulo 

Pres Prudente 

83716 


Bauru 

83722 


Jau 

83723 


Tiete 

83777 


Ataliba Leonel 

83803 


Tatui 

83816 





Table 1. Meteorological Stations Used to Derive Data Sets 
for the Brazil Wheat Model, 





Figure 2. Location of Meteorological Stations Used to Derive 
Data Sets for the Brazil Wheat Model. 


The following Is the selected model: 


OUM 

Dummy variable to adjust Santa Catarina's yield 

TM7 

July Mean Temperature 

P-PET8 

August P-PET 

SP-PET8 

Squared August P-PET 

P-PET9 

September P-PET v 

RDFN11 

Deviation from normal of November precipitation 


Too high temperatures In July and too much rainfall In November both reduce 
yield. The negative coefficient for the linear term P-PET8 suggests that In 
August excess precipitation above demand, PET8, Is damaging to yield. This Is 
reasonable in Brazil when In August the crop is In the tillering stage. The 
quadratic term Indicates that Increased yield Is favorable at some level when 
PET Is higher than precipitation. However, one Is cautioned not to extend this 
Interpretation beyond the limits of the data base used. The statistics of the 
selected model are summarized In Table 2. 

The same variables were used In regression equations for only the two states 
of Rio Grange do Sul and Santa Catarina. The problem encountered was that models 
with agronomlcally reasonable variables had too low an R?; models with 
acceptable R? had too many variables for the number of degrees of freedom. 

Finally, modeling was attempted for the northern states of Mato Grosso and 
Sao Paulo. Overlapping plots of weather variables versus yield were made for 
both states. From these plots. It was determined that modeling for a 
combination of data for these two states would not be acceptable; their climates 
are too different. Next, reasonable weather variables were tried in regression 
for both states separately. No suitable models were derived. For information, 
plots of yields for both states are shown In Figure 3. 


TEST RESULTS 


A jackknife test was run on the selected model. In this test, a year was 
eliminated from the crop data and the model was used to predict that year's 
yield. This process was done for each successive year beginning with 19??„ 
The test had to be run separately on each state. The results are prlr; rd on 
Tables 3 through 5 and plotted on Figures 4 through 6. 
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APPENDIX 


Definition of Variables 

P-PET, precipitation minus potential evapotransplration, Is used a measure 
of the amount of moisture available for plant growth. Potential evapotransplra- 
tlon Is determined by the procedure developed by Thornthwalte (1948). It re- 
quires only temperature: 

PET - X ; 

where I ■ heat index, which is the sum of the 12 monthly indices 1 , 



T 8 monthly temperature in °C, and 

a 8 an empirical exponent 6.75 x 10"^I^ - 7.71 x 10 "^i 2 + 1.79 x io _ 2i + 
0.49. 


The duration of daylight is used to adjust potential e apotransplratlon as a 
portion of 12 hours. 


Statistical Analysis System Procedures Used 

PROC CORR Computes correlation coefficients between 

variables, including Pearson product-moment 
and weighted product-moment correlation. 

PROC PLOT Graphs one variable against another, producing 

a printer plot. 

PROC STEPWISE Provides five methods for stepwise regression. 

Stepwise is useful when selecting variables to 
be Included In a regression model from a 
collection of Independent variables. 


PROC STEPWISE FORWARD Begins by finding the one-variable model that 

produces the highest R*. For each of the 
other Independent variables, FORWARD 
calculates F-statl sties reflecting the con- 
tribution to the model If the variable were 
to be Included. 
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PROC STEPWISE BACKWARD 


PROC STEPWISE STEPWISE 


PROC STEPWISE MAXR 


PROC PETM 
PROC ZINDEX 


Begins by calculating statistics for a model 
Including all the Indepedent variables. The 
variables are deleted from the model one by 
one until all the remaining variables produce 
F-statlstlcs significant at the .10 level. 

The stepwise method Is a modification of the 
forward selection technique, differing In 
that variables already In the model do not 
necessarily stay there. After a variable Is 
added (as In the forward selection method) the 
stepwise method looks at all the variables 
already Included 1i the model and deletes any 
varlab •» that does not produce an F-statlstlc 
significant at the .10 level. Only after this 
check Is made and the necessary deletions 
accomplished can another variable be added to 
the model. 

(Maximum R? Improvement) Unlike the three 
techniques above, this method does not settle 
on a single method. Instead It looks for the 
"best" two-variable model, the "best” three 
variable model, and so forth. 

Uses latitude and mean monthly temperature to 
calculate Thornthwalte's potential evepo- 
transplratlon for each month. 

Uses monthly PET’s, precipitation, SS 
(beginning moisture In surface layer), AWCS 
(available water capacity In surface layer), 

SU (beginning moisture In the underlying 
layer), and AWCU (available water capacity 
In the underlying layer) to calculate Palmer's 
soil moisture budget, drought Index Z, ET, and 
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Figure 4. State of Parana Wheat Model 
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